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BOW WAVE - DETAILS OF ANALYSIS 

By Walter G. Vincenti and 
Cleo B. Wagoner 

SUMMARY 


Details are given of calculations of the flow with detached bow 
wave past a doubly symmetrical, double-wedge profile at zero angle of 
attack. The results of these calculations have been described pre- 
viously in NACA TN 2339* The mixed flow over the front of the profile 
is determined by the solution of a boundary- value problem for the 
transonic small -disturbance eq.uation in the hodograph plane. The solu- 
tion is obtained by the use of finite-difference eq.uations and relax- 
ation techniques. The methods follow established lines except for the 
somewhat novel treatment of the boundary conditions along the shock 
polar and sonic line. The purely supersonic flow over the rear of the 
profile is determined by means of the method of characteristics special- 
ized to the transonic small-disturbance theory. 


INTRODUCTION 


The present report is the second of two NACA publications on a 
theoretical study of the transonic flow past a wedge profile with 
detached bow wave. The specific problem under consideration is that 
of a thin, finite wedge at zero angle of attack, with particular appli- 
cation to the case of the doubly symmetrical, double-wedge profile. 

The first report (reference 1) contains a nonmathematical outline of 
the general analytical method and a discussion of the final calculated 
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?results. The present report supplements the previous work by supplying 
the details of the analysis.^ 

The plan of the report is briefly as follows: In the first section, 

the basic problem of the finite wedge with detached bow wave is stated 
as a boundary-value problem for the transonic small -disturbance equation 
in the hodograph plane. Except for the introduction of a boundary con- 
dition along the sonic line to replace conditions previously prescribed 
in the supersonic portion of the hodograph, this material follows th^ 
lines established by Frankl (reference k) and Guderley (reference 5) • 

It is recounted here primarily for the sake of completeness . In the 
second section of the report, the boundary-value problem in the hodo- 
graph is reduced to a system of finite-difference equations the solution 
of which is then obtained by relaxation techniques. This portion of the 
work, which constitutes the main contribution of the present report, is 
discussed in some detail, since it is anticipated that the methods and 
equations which are presented will be useful in the solution of other 
problems involving detached shock waves. The third section of the 
report describes the transformation of the hodograph solution for the 
finite wedge back into the physical plane. The fourth section is con- 
cerned with the characteristics construction used to obtain the purely 
supersonic flow over the rear of the double-wedge profile, and the final 
section contains a few remarks on the accuracy of the solution. It 
will be presumed that the reader is familiar with the general discussion 
of the investigation given in reference 1 or with the description of the 
finite-wedge problem available in reference 5 , 6 , or 7 .^ 


lAs explained in the introduction to reference 1 , it was originally 
intended that a third report be included under the present general 
title. This report, written by Arthur E. Bryson, Jr., was to describe 
an experimental study of wedge profiles conducted at the California 
Institute of Technology. Because of a broadening of the experimental 
work to include biconvex profiles and subsonic as well as supersonic 
flight speeds , it now seems advisable to depart from the original 
plan. Bryson’s report is therefore being issued by the NACA under a 
less restrictive title (reference 2 ) . Certain of the results of both 
the theoretical and experimental studies have been given in prelimini- 
nary form by Liepmann and Bryson in reference 3 » 

^The authors are indebted to William A. Mersman of the Ames Laboratory 
for suggestions leading to certain of the mathematical procedures used 
in the analysis. 
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NOTATION 



Primary Symbols 

9 . 

critical speed (i.e., speed at which the speed of flow and 
the speed of sound are equal) 

C 

airfoil chord 

e,f ,h,"k 

length of irregular lattice intervals relative to that of 
the basic interval 

F(k),G(X) 

fimctlons defined along sonic line in holograph 
(See equation (6).) 

Iv 

integral defined by equation (39) 

Ji 

(i = 1, 2 , 3) component integrals 
(See equation (33) et seq.) 

K,L 

functions of e, f, and tjq 
(See equation (27).) 


constants 

(See equations (7) and (15).) 

M 

Mach number 

Sin) 

function defined by equation (2k) 

t 

airfoil thickness 

V 

local speed of flow 


Cartesian coordinares 

/v» 

Y 

ordinate finaction 
( See page 31 . ) 

P 

absolute value of at left-hand limit of lattice 

7 

ratio of specific heats (l.k for air) 

A 

basic lattice interval 

n 

holograph variable defined by equation (l) 

e 

local inclination of flow 


k 
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0 ' 


0 


w 


^ /V# 

©u^Qd 


X 


X' 


to 

P 


o 


variable of integration 
(See equation ( 8 ).) 

semiangle of wedge 

ordinates of upgoing and downgoing characteristics at *t) = 0 
( See page 37 . ) 

hodograph variable 
(See equation ( 6 ).) 

variable of integration 
( See equation ( 6 ) . ) 

speed parameter 

(See equation ( 49 )») 

transonic similarity parameter 
(See equation (l 7 ).) 

fluid density 

stream function 

value of ^ at the point E 
(See equation (16).) 


Subscripts 


conditions in free stream 


o,ija,etc. value at a prescribed lattice point 
* conditions at critical speed 


Superscript 


quantity in normalized form 
( See equation ( 9) • ) 
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STATEMENT OF BOUNDARY-VALUE PROBLEM IN 
HODOGRAPH PLANE 

Complete Hodograph 


A qualitative representation of the flow over a double-wedge 
profile with detached bow wave is given in the first of the accompanying 
sketches, which differs only in detail from a similar sketch presented 
in reference 1. A description of this flow field is to be found in the 



Shock waves 

Sonic tine (V/a,= !) 

— Streamline 

Expansion \ Mach 

Compression ) lines 


earlier paper. The second sketch shows the corresponding picture for 
the flow about the front half of the profile in the hodograph plane - 



that is, in the plane V/a*, 6 where V is the local speed of flow, 
a* is the critical speed, and 6 is the local inclination of flow 
relative to the x axis. 


The picture in the upper half of the hodograph plane can be 
described briefly as follows; The part of the shock wave which borders 
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upon the subsonic region in the physical plane appears in the holograph 
as the subsonic portion AE of a shock polar. The shape and position 
of the shock polar are determined by the dimensionless free-stream 
velocity (or, what is equivalent, by the free-stream Mach number 

Mq) and by the ratio of specific heats 7.® The portion of the central 
streamline from the normal part of the shock wave to the stagnation 
point at the nose of the wedge maps into the portion AO of the hori- 
zontal axis in the holograph. The image of the wedge itself is given 
by a radial line inclined at the wedge angle 0-^ and extending from 
the origin 0 to the point B on the critical circle (v/a* = l) . 

The shoulder of the wedge, which produces an expansion fan of a locally 
Prandtl-Meyer type in the physical plane, appears in the holograph as a 
portion of the downgoing characteristic (epicycloid) starting at B. 

The last Mach line from the shoulder to the sonic line (termed the 
separating Mach line in reference l) appears as part of the upgoing 
characteristic which begins at the intersection E of the shock polar 
and critical circle. Point G, the point of intersection of the epi- 
cycloids from B and E, fixes the extent of the downgoing character- 
istic which must be considered in determining the solution in the 
holograph. A typical streamline in the holograph plane is shown by the 
line FCD. 

To obtain a solution for the detached-wave problem in the holograph, 
a boundary -value problem for the differential equations of gas dynamics 
must be solved within the region AOBGEA. If the stream function ^ 
is taken as the unknown, the pertinent boimdary conditions are as 
follows: 

1. The value of >|r is constant along the basic stre aml ine AOBG. 

2. The streamlines (i.e., the lines of constant leave the 
shock polar with a direction which is a known function of location on 
the polar. 

3. The increment in ^ over the portion AE of the shock polar 
has a prescribed value different from zero. 

The reason for the first condition is obvious. The second condition is 
a consequence of the requirement that, at every point on the shock polar, 
the direction of the shock wave as computed from the solution for ijr 
must be compatible with the direction given by the equations, for an 
oblique shock wave in a uniform stream. The third condition prevents ^ 
from being simply a constant throughout the holograph and, in effect, 
fixes the scale of the flow field in the physical plane. It will be 
noted that no condition is prescribed along the boundary EG in the 
holograph. Frankl has proved (reference k) that the solution determined 
by the foregoing boundary conditions is unique. 


®The equations which are pertinent here can be found in the work of 
Frankl (reference . 
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Specialization to Small Disturbances 


Original boundary-value problem .- As has been shown by Guderley 
(reference 5) > the equations of the boundary -value problem in the holo- 
graph are considerably simplified when restriction is made to the 
neighborhood of the critical speed. To this end, the quantity is 
introduced according to the relation 


r, = (,+i)V3 

U 5|C 


( 1 ) 


and T) and the stream angle 0 are assumed sufficiently small that 
only their lowest powers need be retained in the analysis. This means, 
in effect, that the right-hand portion of the previous holograph 
(including the shock polar itself) is made to shrink down to the vicin- 
ity of the point H, which defines the intersection of the critical 
circle and the horizontal axis. 

When the foregoing procedure is carried out and the limiting 
process is counteracted by a suitable enlargement of scale, the situ- 
ation in the small -disturbance holograph (i.e., in the t], 0 plane) 
appears as in the adjacent sketch. Here the critical speed corresponds 





to the vertical axis t] = 0. The equation for the upper half of the 
shock polar in the simplified holograph has the form 





(2) 
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where t]q is the value of t) corresponding to the free-stream 
velocity Vq.^ By virtue of the limiting process, the stagnation point 
at 0 has moved, in the present system of axes, infinitely far to the 
left. As a result, the part AO of the horizontal axis ($ = O) extends 
nqw from t] = -tiq to t] = The image OB of the wedge is similarly 

represented by the horizontal line .0 = 6-^, t] ^ 0. The characteristics, 
which complete the boundaries of the field, have the simple form 


6 = const. 


± - -n®/^ 
3 


( 3 ) 


On the basis of the usual assumptions regarding flow near the 
critical speed, the differential equation for .\(r reduces, in the pres- 
ent simplified hodograph, to the form 


^TJT) ■ 

This is the mixed elliptic -hyperbolic equation studied by Tricomi in 
reference 8. The boundary conditions along the central streamline 
require that ^ is constant - say 0 - on AO and OBG and that 

\|r->0 as T]->( -oo) for 0 ^ 0 ^ 0w 'the upper half of the shock polar, 

the boundary conditions require that the lines of constant t must have 
the slope 


^ = _ /no+Tl 

dq 3 t1q+5i1 ^ 2 


(5a) 


On a line of constant >1^, d0/dq can be replaced by so that 

the fofegoing condition can also be written 


. jtoiTn . 0 (5b) 

The final boundary condition requires that ''If must have some given 
value at the point E. Since the coordinates of the flow field 

will ultimately be e3q)ressed in terms of a characteristic dimension of 
the wedge, the actual value assigned to is purely a matter of con- 

venience. As before, no boundary condition is specified along the 
characteristic EG. 


'‘'The derivation of this and the other equations for the simplified 
hodograph is given by Guderley in reference 5* 
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Elimination of the supersonic region .- The foregoing is the 
boundary-value problem for the finite wedge as formulated by Guderley. 

It was the original intention in the present work to obtain a numerical 
solution of this problem on the basis of the boundaries and boundary 
conditions which have been described. Efforts in this direction failed, 
however, because of difficulties in obtaining convergence of the relax- 
ation process in the supersonic portion BGE of the holograph.^ 

Similar difficulties have been reported in references 9 and 10 with 
regard to relaxation calculations of the transonic flow through a 
converging-diverging nozzle. The reasons for the difficulty in the 
present case are not apparent. Fundamental questions would appear to 
be involved concerning the stability and convergence of the finite- 
difference scheme for the Tricomi equation in the hyperbolic domain. A 
study of these matters, similar perhaps to that reported for the wave 
equation in reference 11, may be a prerequisite to numerical solutions 
of mixed-flow problems in the general case. In the present example, 
however, the difficulty can be circumvented by modifying the boundary- 
value problem so as to eliminate the supersonic region from explicit 
consideration . 

The elimination of the supersonic region depends on a formula 
given by Tricomi (reference 8, equation (2.19)) which relates the 
behavior of ^ on the vertical axis to its behavior on a character- 
istic. In the present case, in which ^ is identically zero on the 
characteristic BG, this formula reduces to an integral relation 
between ^ and \|t«q at points on the sonic line. This relation has 
the form 


F{X) 


k, fV gjill- 


dX,' 


( 6 ) 


where 


X = Q-\ff “ S 

f(x) = tCo^e) 

G(x) = tT^(o,e) 

and X' is a variable of 
integration. (See sketch.) 





^Several procedures were tried in the supersonic region, using both a 
square lattice and a lattice following the characteristics. All were 
unsuccessful. 
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The niimerical constant is given by 


^ 3^/3 r°(i/ 3 ) (T,) 

4jt2 

where r(l/3) is the gamma function of the argument l/3« For the 
accuracy required in the later numerical work, the value of this con- 
stant may be taken as 

kj^ = 1.013 (Tb) 

Satisfying equation (6) everywhere on the sonic line from B to E is 
completely equivalent to satisfying the condition i = 0 on the down- 
going characteristic from B to G. 

For the present application, it is convenient to invert equa- 
tion (6), which can be done by means of Abel's formula. This gives 
(see reference 12, p. 229) 

g(m = -ZIa r" 

2nkj^ dX ^ 

The differentiation indicated on the right is readily accomplished by 
first transforming the integral to one with fixed limits by means of 
the substitution X' = tX. The result is, after reverting to the orig- 
inal notation. 


G(X) = - 


Transforming 
that F(o) = 


JT 1 

" 1 


r’' a>,l 

2 itki \ 

_ 3 

J (X-X')2/3 

J (X-X')2/3 

o ' 


the first integral through integration by parts and noting 
0, one obtains finally 


G(X) 
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This can be written in the r|,0 notation as * 



d.0' = 0 


(8) 


where 0*= denotes the variable of integration. As with 

equation (6), satisfaction of equation (8) everywhere on the sonic 
line between B and E insirres that ^ is zero everywhere on the 
characteristic from B to G. By regarding equation (8) as a 
boundary condition along BE, the region of solution of the partial 
differential equation (4) can be confined to the purely subsonic 
portion of the holograph (11*^0). Relaxation methods can be used to 
solve the resulting elliptic problem without essential difficulty. 

Equations in normalized form. - To carry out the numerical 
calculations, it is convenient to normalize the equations of the 
boundary-value problem by means of the transformation 


This is equivalent to introducing the rules for transonic similarity. 
(See, for example, references 5^ 13^ and 1^.) The particular form 
of transformation chosen here has the advantage for the present work 
of providing a unique shock polar with conveniently located hori- 
zontal and vertical intercepts. 

With the foregoing substitution, the differential equation (4) 
takes on the following form in the plane: 



e 


(9) 



(10) 


\ 


L 
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Consistent with the elimination of the 
value problem can now be summarized as 

supersonic region, the boundary, 
follows ( see sketch) ; 



§ 



B 








1. On the basic streamline AOB: 


t = 0 for 0 = 0 , 't] = -1 

(11) 

ik = 0 for 0 = 0^, r)^0 

(12) 

^ 0 for T) -»-oo^ 0'^ 6 ‘^6-^ 

( 13 ) 

On the shock polar AE: 


‘/l+n - 0 

T] 3+'5T1 0 



for 


0 = (i-'Ti)yi+?f, - 


m 


MCA TN 2588 


13 


3. On the sonic line BE: 


+ kg 

T 1 



a ♦g(o,e') 


for 



where 



(15) 


At the point E: 

= ikg (arbitrary constant ^ 0 ) for ‘1^ = 0, "^=1 


(16) 


It is apparent from the preceding equations that a solution of the 
problem will depend on only the single p^ameter which defines the 

position of the upper boiindary in the r\,Q plane. This parameter is 
directly related to the more familiar transonic similarity parameter 


which was used as a basis for the presentation of the results in 
reference 1 . (in this equation, t/c is the thickness ratio of the 
double-wedge profile. It is equal, to the present order of accuracy, 
to the wedge angle G^.) The relation between and 0 ^, easily 

derived from equations (l) and ( 9 ) of the present paper and equa- 
tion (21) of reference 1 , is given by 


The solution of the boundary-value problem in the t],0 plane is 
obtained in two steps, according to established procedures for the 
numerical treatment of partial differential equations. (For intro- 
ductory articles, see references I5, 16, and I7. For an extended 
discussion, see reference I8.) In the first step, the domain under 


[(7+l)(t/c)]^/® 



(17) 


21/3 
O ~ -^^2/3 


(18) 


SOLUTION OF BOUNDAEY-VALUE PROBLEM 
IN HODOGRAFH PLANE 
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consideration is covered by a square lattice, and a finite-difference 
approximation to the differential equation or boundary condition is 
vritten for each lattice point. The boundary -value problem for the par- 
tial differential equation is thus reduced to a problem in solving a 
large number of simultaneous algebraic equations. Solution of the 
latter problem by relaxation methods is the second step. 


Reduction to Finite-Difference Equations 


The arrangement of a typical finite-difference lattice in the t\,G 
plane is shown in figure 1. The basic lattice interval, which is the 
same in both directions, is denoted by A. Adjacent to the shock polar, 
the interval is adjusted so that the terminal lattice points lie on the 
polar itself. For purposes of formulating the finite-difference equa- 
tions, the lattice points are conveniently grouped into five categories 
as follows (typical points in each category are indicated in the figure) 

a. Regular internal points 

b . Points far to the left 

c . Points adjacent to the shock polar 

d. Points on the shock polar 

e. Points on the sonic line 

The form of the finite-difference equation pertinent to each 
category will be developed in the subsequent paragraphs. The methods 
employed are standard, except for the somewhat novel treatment of the 
boundary conditions along the shock polar and sonic line. 


Regular internal points .- The category of regular internal points 
conqprises all points interior to the boundaries but not immediately 
adjacent to the shock polar. The situation in the vicinity of such a 

point is as shown in the sketch. The differ- 
ence equation which applies here is obtained 
by suitable approximation to the differential 
equation (lO). 


4 



A 


If it is assumed that the unknown 
function ^ niay be expanded locally 

in the form of a Taylor’s series, the values 
of \lr at 2 and 4 may be written (see sketch) 
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^2 


1^0 + A 


o 


+ — 

21 


T)T) 


O 



T]T\T] 


+ o(a*) 


^4 = to - At^l 
T]|o 



^ I 

Tn|o 


31 


+ 0(A*) 

TlTlTljo 


Addition of these equations and solution for gives 

n n|o 


i — 1 

TIT) |o 


- (t4-2to+t2) + 0 (a2) 
A^ 


(19a) 


which is a well-known difference exj)ression for the second derivative. 
The corresponding derivative in the vertical direction is similarly 
represented hy 


= h ('|'3-2v^i) + ( 19 b) 

6 9 o ^ 

Substituting these expressions into the differential equation (lO) and 
neglecting the terms of O(a^) then gives for the finite -difference 
equation at a regular internal point 

>1/2 + - 2T)o(t 1+11^3) - 2(l-2?)o)i^o = 0 (20) 

where rfo denotes the horizontal coordinate of the point in question. 
The difference equation for internal points is thus the same for points 
on a given column but differs from one column to the next. For a point 
adjacent to the upper boundary OB — as, for example, the point a* in 
figure 1 — the value of must be set equal to zero in accord with 

the boundary condition (l2). Similar considerations hold for points 
adjacent to the lower boundary OA. 

Equation (20) represents the simplest possible finite-difference 
approximation to the differential equation (lO). As is apparent from 
the derivation, the error involved is of O(A^) . Consideration has 
been given to improving the approximation by including additional lat- 
tice points in the finite-difference equation or by incorporating 
higher-order difference corrections in the later relaxation work (cf., 
references I9, 20,. and 2l) . Because of the complicated nature of the 
boundary conditions along the shock polar and sonic line, however, 
consistent application of these procedures did not appear feasible. 

The requisite accirracy in the present work has therefore been achieved 
by suitable decrease in the mesh interval A in those regions in which 
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the function ^ varies most rapidly. This procedure has the secondary 
advantage of providing closely spaced values of the derivatives which 
are required for the later transformation to the physical plane. 

Points far to the left .- In order to carry through the numerical 
analysis^ it is necessary that the finite-difference lattice be termi- 
nated at some distance to the left in the holograph. This can be done 
with the aid of an asymptotic solution valid for large negative values 
of rf. 

By separation of variables, it can be shown that the general 
solution of the differential equation (lO) in the region ^^-1, 

0=0= 0^^, subject to the boundary conditions (ll), ( 12 ), and ( 13 )^ is 


^ ^ An sin K ^/3 

n=i 

where 73 is the modified Bessel function of the second kind of order 
1/3 (notation of reference 22) and An is an appropriate constant. At 
sufficiently large negative values of q the first term of the series 
will predominate, and the above solution can be approximated by 


^ (- 2 ^) 3/2 
L3ew 


^ = A 


sin 



4- 

L30W 


If ^ 1/3 is then replaced by the first term of its asymptotic expansion 
(reference 22, p. 202) 



there results finally for '|f the expression 


t = B sin ^ ^ 

) X ( -^)"^/'^ exp 

4- (-2'^i)^/^' 



L 30W J 


( 21 ) 


where B is an unknown constant. 

The asymptotic solution (2l) makes it possible to terminate the 
finite-difference lattice at a position on the left. Consider a typi- 
cal lattice point in a column located at rJ=-0 (as, for example, the 


HA.CA TN 2588 


17 


point b in. fig. l) . The neighboring points 
are then as shown in the adjacent sketch, where 
the point k now represents a ficticious lattice 
point located at tT = - ( 3+A) . If 3 is 
taken sufficiently large that a/P<< 1, then 
it follows from equation (2l) that, to a first 
order of approximation, ^ • 



Substitution of this value of into the 

previous equation (20) gives for the finite -difference quation at a 
point on the left-hand boundary 



+ 23(\|f^+t3) + exp ^ - 2(l-h23) ilfg = 


( 22 ) 


The value of 3 to be used in any particular case is determined on the 
basis of computational experience. In general, the larger the value 
of "the larger must be the value of 3 "to assure that the use of 

the asymptotic solution (21) is justified. Since the over -all result 
is insensitive to changes in the left-hand portion of the field, how- 
ever, the choice of 3 is not a critical matter. 


Points adjacent to the shock polar.- Points adjacent to the shock 
polar require special treatment because of the irregularity of the 
intervals encountered near the curved boimdary. Consider the typical 
case shown in the sketch (corresponding to point c in fig. l) . 

Here h and k define the length of the 
irregular intervals relative to that of the 
regular interval A. 


To obtain the desired accuracy, it was 
found advisable in the present case to 
include three rather than two neighboring 
points in each of the coordinate directions. 

The value of t at the points 2, k, and 6 
is therefore written 

6 


/ 
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Ik = \|f + kA 'tr I + ■ t 

2 O ^ P' 01 

T1 0 “=• TIT p J* 11 11 T) 


(kA)® 

+ — \|r 


e\j r\^ 


+ O(A^) 


1114 = ikn-A to 


+ — 


0 2i T] T) 

0 3 . TITITI 


+ O(A^) 


( 2A) ^ 

tg = to - 2A i + ^ t 


( 2 A)' 


^ r\ • -w-^l • *-w I 

r\\o 21 Tj T] io 3* T1 T |o 


^ 1 + O(a^) 


These may be looked upon as constituting three simultaneous equations 
for the first three derivatives of t in the horizontal direction at 
the point 0. Solution of these equations for t,v,~ gives 

T)T) - 


t, 


■ ^ - - *s * 

r\Vi\o L 2+k 


t. - — tn + t 

1+k ^ k ° k(l+k)(2+k) 


+ o(a^) 


The corresponding expression for t~~| is identical except that k is 

6 010 

replaced by h and tg , t4, and % by tg, ti, and ts, respectively. 
Substituting these expressions for the two second derivatives into the 
differential equation (lO) and neglecting terms of O(a^) then gives 
for the finite-difference equation at 0 


6 2(2-k) . 1-k , 

k(l+k)(2+k) ^ i+k ^ ~ 2^k 


6 


h( 1+h) ( 2+h) 


^3 + 


2(2-h) 

1+h 


\lTi - ijf ! 

^ 2+h 



^0=0 


(23) 


This reduces to the previous equation (20) when k = h = 1. (The 
functions of h and k which appear here have been tabulated in 
reference 23. The intervals of tabulation are not always sufficiently 
small, however, to provide the accuracy needed in the present work.) 

Points on the shock polar. - In past applications of numerical 
methods to problems involving curved boundaries, it has not ordinarily 
been the practice to use a lattice with points located on the boundary 
itself. The prescribed boiindary conditions have then been incorporated 
in the following manner (cf. references 2 h and 25): First, the finite- 

difference lattice is extended, on the basis of the regular lattice 
spacing, to include ficticious points external to the boundary. This 
makes the lattice geometry at internal points adjacent to the boundary 
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the same as at all regular internal points. Next, with the aid of the 
boundary conditions and suitable interpolation and extrapolation formu- 
las, an expression is obtained for the independent variable at each 
external point in terms of the values at neighboring internal points. 
Finally, by substituting these expressions into the finite-difference 
equation for a regular point, the difference equations are written for 
the internal points adjacent to the boxmdary. In this way the boundary 
conditions are incorporated implicitly into the difference equations at 
internal points. The procedure is parallel in many respects to that 
used in terminating the present lattice at the left-hand side of the 
field. 

Although a procedure of the foregoing type can be devised to take 
care of the boundary conditions on the shock polar, a different approach 
was found advantageous for the present work. In this approach, the 
lattice points are placed directly on the boundary as previously 
described, and a difference equation is obtained at each such point by 
suitable finite-difference approximation to the boundary condition. 

This leads to a somewhat larger system of simultaneous equations than 
would the more usual procedure but has been found in the long run to 
give more accurate results with less total effort. 

The boundary condition (l4), which is thus the basis for the 
finite-difference equations on the shock polar, can be written 


- S(^) = 0 (24) 

11 e 


where 



The problem now is to determine differ- 
ence expressions for the derivatives 

and at points on the polar, 

r) e 

Consider, for example, the typical 
situation shown in the sketch (corre- 
sponding to point d of fig. l) . To 

determine the value of t at each 

0 

of the interior points 1, 2, and 3 is 
written, as before, in terms of a 
Taylor’s series about the boundary 
point 0. Solving the resulting three 
equations for the derivative at 0 

then gives d 
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t. 


2 ^ 

A 


_ h(2+h)+2(l+h)^ ^ ^ (l+h)(2+h) 

h( 1+h) ( 2+h) ° 2h 


1+h 2( 2+h) 


+ o(a®) 


(25) 


This expression, which includes terms of O(a^) > is inconsistent in 
order of acciaracy with the expressions previously employed in setting 
up the finite-difference equations at internal points. Since end 
differentiation is, even for a given order of mathematical accuracy, 
inherently less precise than differentiation at a midpoint (see error 
terms in reference 23) , the retention of the second-order terms was here 
thought advisable. 

The determination of the corresponding expression for is a 

bit more involved. Expanding ilf at the boundary points U and 5 
means of a Taylor's series in two dimensions gives 


4-4 = ^0 - A i~l 

- fA t~| 

A^ 
+ -~- 


+ 


1 + 

(fA)2 

t<v.v I 

+ 

nl 

|o 

jo 

21 

n^i! 

0 

Tie 

|o 

21 

0 0 1 

0 

tg = to + A t,. 1 

1 + eA 1 

A2 

+ 


1 + 

eA^t 1 

1 + 

CM 

^ c 

t'. ^ 1 

1 + 

Tl' 

|o 

G 10 

21 

11 T) 

|o 

qe 

|o 


00 

lo 


An expression for is already known from equation (25) in terms 

e )o 

of ■'1^0 j ^35 expression for can similarly be 

© 0 lo 

determined. The two foregoing expansions may thus be regarded as con- 
stituting two equations for the three unknowns • 

n'o niilo nelo 

To solve for t^l , one more equation is necessary. This is provided 

T)|o 

by the differental equation (lO), which also applies on the boundary 
and which may be written at the point 0 as 


- 2qo'^-^.-l =0 

Tin o eelo 


The solution for is then found as 

qjo 
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j _ 

1- ~ K 

_ l|fQ + .JL \|ft. 

- ^ t 1 

jo A 

t e+f 

e+f e+f 

e+f e 1^ 


1 e-f 


§_ to - ^ t. 

h(2+h) ° h ^ ^ 2+h '• 


O(A^) 


(26) 


vhere ^»w| is given by the previous equation (25).® 

e |o 

The required finite -difference equation for the point on the 
shock polar can now be obtained by substituting expressions (25) and 
(26) into the boundary condition (2k) and neglecting the higher -order 
terms in each case. The result can finally be written 


(l+h)(2+h)+2L(3+2h) , h(2+h)+UL(l+h) 

t. 


2h 


1+h 


h(l+h)+2L(l+2h) 

2(2+h) 


*3 + 


eK - fK tg - 


h( 2+h)+2(l+h)^ 
h( 1+h) (2+h) 


6l 


h( 2+h) 


+ (e-f)K 


to = 0 (27) 


where K = K (e,f,T)o) and L = L (e,f,qo) are given by 


(e+f) 3(710) +2ef 


L = (e-f) (^Tio - K 


Equation (27) is convenient for points on the shock polar for 
which -0.6 ^ T)<0. For the general procedure is the 


®It will be noted that the coefficients in equation (26) become 

undesirably large as f — >-e and are undefined when f = -e. This 
results from the fact that the determinant of the coefficients in the 
simultaneous equations used to obtain vanishes when f = -e. 

Difficulties from this source can be avoided by judicious choice of 
the lattice points. 
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same except that the points 
0, 1, 2, and 3 ai’e now more 
conveniently located on a hori- 
zontal line and the quantities 
e, f, and h are redefined 
accordingly ( see sketch) . The 
resulting finite -difference 
equation is identical with 
equation (2?) except that the 
terms which previously arose 

from the expression for t*-* 

0 

are now multiplied by -1. 

Points on the sonic line .- 
The difference equation for 

points on the sonic line is obtained by finite-difference approximation 
to the boi^dary condition (15), which can be written, to the accuracy 
required in the numerical work, as 



-f 0.342 / ® 77 ^- .c --r 7 

T) J 0' -0)2/3 


d0' = 0 


(28) 


The procedure varies depending on whether the difference equation is 
being written for the first point below the upper boundary or for one 
of the lower points. 


First point below boundary: 

At the first point below the upper boundary the situation is as 
indicated. To approximate the integral in equation (28), use is made 

of Guderley' s singular solution for tran- 
sonic flow over a convex corner (refer- 
ence 26) . On the basis of this solution, 
it can be shown that the variation of ilr 
along the vertical axis in the immediate 
vicinity of the point B is of the form 


B 


// 


O 




(29) 


If the lattice spacing is made sufficiently small, this asymptotic 
relation may be taken as approximately correct over the entire interval 
from B to 0, so that within this interval 
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and 




/^y/3 

A J 


9 


4 (9^-9)^/^ 


For "the first point beiow the boundary, 
can thus be written 


the Integral in equation (28) 



t§r(o,g'’) 

(0- -0)2/3 


d0' 


***0 rS^-A 

J, [A-(V8’ 

®w 



t1/3 


dr 


where t = (0v-^»)/a. 
standard form by means 
gives 


The integral on the right can be reduced to 
of the substitution t(1-t) = i/4 which 


^ dT = r " 

(1-t) 3/3 2^/3 y 

This is an elliptic integral of the first kind. Its value, as deter 
mined from the equations and tables of reference 27, is 



o 


p 1 r ^/3 

J U-t)3/3 


3 ^ 1.846 

22/3 ^1/4 


2.650 


The integral in equation (28) thus becomes, in the present case 



>|f^(O,0') 

(“S'-?) 


d0' 


4 

- X 2.650 

3 


^ o 
^2/3 


3.533 


lo_ 

2,2/3 


(30) 
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To approximate the derivative in equation (28), 4' is 

T1 

expanded at points i and ii by means of a Taylor's series about 
point O. Terms involving may be omitted here, since the 

Ti Tito 

differential equation (lO) shows this derivative to be identically zero 
at points on the sonic line. The values of at i and ii can 
thus be written 


- A 


— ^ \1/’ 

o 3'. mn 


+ o(a^) 


I|r. . = T|r^_ 2A t~| - — 1 + O(A^) 

° Tilo 31 n mlo 


Solution of these equations for gives, to the second order in A, 

T] O 




(31) 


Substitution of expressions (30) and (31) into the boundary condi- 
tion (28) gives the following finite-difference equation for the first 
point below the upper boundary: 


A + 1.208 A^/ 3 ^ ^0 = 0 (32) 

3 6 u 

Unlike the previous equations (20), (22), 
involves the value of A. 


(23), 


and (27) this equation 
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Lower points: 


For a general point below the first lattice point, the situation 
is as represented in the sketch. The integral in the boundary condi- 


tion (28) is here evaluated in three sections. The 
integral over the lattice interval from B to n -1 
is evaluated on the basis of the asymptotic relation 
used before. The integral from n -1 to 1 is 
evaluated by assuming a linear variation in \|r over 
each of the included Intervals and then integrating 
analytically. The linear assumption is sufficiently 
accurate here, since this entire middle section con- 
tributes only a relatively small portion of the 
complete integral. The integral from 1 to 0 , 
which contains a singularity in the integrand at 
the point 0, is evaluated by expressing as a 
cubic in terms of its value at the points 2, 1, 0, 
and -1 and then integrating analytically as before. 
The added accxiracy of the cubic is required here, 
since this last section contributes by far the 
majority of the over-all value. The boundary condi- 
tion (28) may thus be written 


( 0 , 9 ) + 0. 3*^-2 (J1+J2+J3) = 0 


where the J's represent the three component 
integrals just described. 


(33) 


B 


n 


n-1 


n-2 


mi-/ 


m 


// 


-/ 


Proceeding to the details of the above procedure, the integral 
from B to n-1 is first written 


Ji = 


t'§( 0 , 9 >) 


^W 


[nA-(V9')]^/= 


d 9 ' = 


4t. 


n -1 


jl /3 


3A‘ 


2/3 




dr 


where, as before, t = ( 9 ^- 9 ' )/A. This integral can be expressed, if 
desired, as the difference of two elliptic integrals of the first kind. 
For present purposes, however, it is more convenient (and sufficiently 
accurate) to expand (n-r)"^/® according to the binomial theorem and 
integrate termwise. This gives finally 


Jl = 


(nA) 2/3 


, 8 2 

1 + + + 

21n 9n2 


160 


55 


1053n® 486 n^ 


t 


n-1 


(34) 


This expression is used, of course, only for n^ 2 . 
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On the basis of the assumption of a linear variation of 't between 
adjacent lattice points, the integral from n-1 to 1 becomes 


Jg 


'0v-(n-i)A 


w 




n-2 


'^m+i" ^m 


m=l 


p6^ -(n-m)A 
0^-(n-m-l)A 


[nA-(0w-0')]^/° 


d0' 


Carrying out tne integration gives 


n-2 



m=l 






(35a) 


This expression is valid for n =3* (For n = 2, Jg obviously does 
not exist.) For n > 3 it is convenient to rewrite the summation so 
that the value of at a given point is not repeated in successive 
terms of the series. This is done by separating expression (35a) into 
two series (one with and one with i) > expanding these series, 

and then regrouping terms. The result is finally 


^2 = 


. 2/3 


2i/a-i 


n-2 


m=2 


(m+l) ^/®-2m^/^ + (m-l)^/® 




(n-l)^/°.(n-2) 





n-1 


} 


(35b) 


This expression is valid for n ^ 

To evaluate the integral from 1 to 0, t is represented within 
this interval by a cubic of the form 
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where a, b, c, and d are determined such that i|r has the proper 
values at the points 2 , 1, 0, and -1. This expression is to be substi- 
tuted into the integral 


J 


3 


6-^-nA 


©w -(n-i)A 




de» 


The result is finally, after evaluation of the coefficients a, b, c, 
and d, 


J 


3 


-1 



_ §z ^ + 33 ^ 

28 14 



(36) 


The finite -difference equation for a general point on the sonic 
line can now be obtained by replacing the J^s in equation (33) by the 
expressions (3^)^ (35)^ (36) and using the previous expression (3l) 

for the derivative . The result is a lengthy linear equation 

T ) I o 

involving the value of \|r at the points i, ii, -1, 0, 1, 2, 
n-1, n. Fortunately for the later relaxation work, the coefficients 
of the terms turn out to be relatively small for all points above the 
point 2. 


Distribution of mesh points.- When an attempt is made to solve the 
present problem with a coarse mesh, it is soon found that most of the 
variation in t takes place in a relatively small region near the 
intersection of the shock polar and sonic line. To obtain a suffi- 
ciently accurate solution in a practicable length of time, it is 
therefore necessary to employ a graded lattice, that is, a lattice 
which has different spacing in different parts of the field. Figure 2 
shows the distribution of lattice spacing found satisfactory in a 
typical case (0^ = 1.6). The particular arrangement shown here 
involves a total of 228 lattice points. For other values of 6 ^, the 
grading of the lattice follows the same general scheme. Obviously, 
however, the total number of lattice points must be increased as the 
upper boundary is moved farther from the shock polar. 

Formulas (34), (35), and (36 ), which are used to approximate the 
integral along the sonic line, presume the existence of lattice points 
at a uniform interval over the full distance from the upper boundary 
to the point in question. This condition is not fulfilled in a graded 


'^Occasionally, when two points on or adjacent to the shock polar fall 
very close together, one of the points is arbitrarily omitted. An 
expression for the omitted value of , which is then necessary to 
complete the difference equation at neighboring points, is found by 
parabolic interpolation between the values at the available locations. 
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lattice such as that indicated in figure 2. Some modification of the 
me'thod must therefore be made to obtain the finite-difference equation 
for a point on the sonic line in one of the regions of finer mesh. 

This requirement was satisfied by means of a simple averaging process 
in which the contribution of nonexistent fine-mesh points is replaced 
by an average contribution expressed in terms of at bracketing 
points on the available coarser net. Since the contribution of indi- 
vidual points is small even for points only moderately removed from 
that at which the equation applies, a rather crude averaging process 
is sufficient in most cases. (The details need not be given here as 
they would soon become apparent to anyone working with the method.) 

When the averaging procedure would not be sufficiently accurate (as 
when the point at which the equation applies is near the line of demar- 
cation between two different sized meshes), fictitious intermediate 
points are introduced into the coarser net and the value of t at 
these points is obtained from plots of the distribution of along the 
sonic line. 


Solution of Finite-Difference Equations 


By the methods of the foregoing section, a finite-difference 
equation can be obtained for each lattice point in the hodograph plane. 
The result is a large number of simultaneous algebraic equations 
involving an equal number of unknown values of \|r. Since the number 
of unknowns in each equation is small, the equations lend themselves 
well to solution by relaxation techniques . ® 

The mechanics of the relaxation process have been well described 
by various authors (references 15, l6, 17, and l8) and need not be gone 
into here. For present purposes it was found satisfactory to take ilfg 
in the boundary condition (l6) equal to 10,000 and work with integer 
values of i|r throughout most of the field. The residuals in the relax- 
ation process were eliminated to within limits of ±2 (with due cafe, 
of course, that all residuals in any given area were not predominately 
of the same sign) . To obtain satisfactory smoothness of the solution 
near the left-hand boiindary in some examples, it was necessary in this 
region to work with values of ijr to 0.1 and eliminate residuals to 
within ±0.5* Whenever the coefficients in the finite -difference 


®It is interesting to note that, of the complete set of simultaneous 
equations, only two — those for the points on the shock polar and 
sonic line immediately adjacent to the point E — are not homo- 
geneous. Only this fact prevents the solution of the complete set 
from being identically zero. 
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equations were relatively small, the corresponding terms were neglected 
in the point-by-point adjustment of il/”. The error so introduced was 
eliminated periodically by recomputing the residuals using all terms in 
the finite-difference equations. This procedure was particularly help- 
ful in the case of the lengthy equations which apply at points on the 
sonic line. The transition between the various regions of the graded 
lattice, which is not often discussed in the literature, was accom- 
plished by the use of overlapping fields in essentially the manner 
described in reference 28. 

By means of the foregoing procedures, the boundary-value problem 
in the holograph plane has been solved for values of \ of 1.3, 1.6, 
2.4, and 4.2. These values correspond, respectively (see equation (18)), 
to values of the transonic similarity parameter of I.O58, O.92I, 
0 . 703 ^ and 0.484 as given previously in reference 1. As shown in the 
earlier paper, calculations for these values were sufficient to bridge 
the gap between the previous results of Guderley and Yoshihara 
at 5o = 0 (reference 29) and the analytic results which can be obtained 
when the bow wave is attached and the flow is completely supersonic 
(reference 30). As an example of the solution in the holograph plane, 
the variation of for = 1.6 is shown as a function of and 
G in figure 3* (These results correspond to the results shown in the 
physical plane in fig. 2 of reference 1.) This figure shows clearly 
the rapid variation of ^ near the intersection of the shock polar 
and sonic line. The calculated values of ^ corresponding to figure 3 
are listed in table I at the end of the report. 


TRANSFORMATION TO PHYSICAL PLANE 
Flow Field 


The transformation from the holograph plane to the physical plane 
is governed, in the small-disturbance theory, by the following equa- 
tions (cf. reference 5): 


dx = 


(7^1) 1/3 


(n te dTi + ^T) de) 


dy = — ^ dll + ^0 d0) 


-i— dt 


The second of these equations implies that, in a flow field determined 
according to the transonic small-disturbance theory, all streamlines 
appear as straight lines parallel to the horizontal axis. When 
expressed in terms of rf and 9 the foregoing equations become 
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The length Z of the wedge, which is equal to one-half the chord 
of the double -wedge profile, can be found by integrating the first of 
these equations over the upper boundary OB of the hodograph. (See 
sketch.) This gives 


e 


o 


a -if. 




I 


c 

2 


(y+1) 

p*% 




*\0 . ^ . tyj 

0 


With this relation the previous transformation equations can be put in 
the dimensionless form 


f - V ^ ( 2rf dT) + d.0) 

^Iw 0 T1 

(7+1)^/® = J:_ d^|r 

V 2 Vc / 


where Iv represents the integral 


Iw = 


~ I / ^ Of \ 

n K dT) 

0 


(37) 

(38) 


o 


(39) 
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To obtain the flow field in the form given in reference 1, equation ( 38) 
must be rewritten in terms of the ordinate function 


Y = [(7+l)(t/c)]^/® (I 


The result, derived with the aid of the relation 
is 


dY ^ [(7+l)(t/c)]^/® " ■— ( 40 a) 

Integrating this relation, subject to the condition that Y = 0 when 
^ = 0, gives 

Y ^ [(r.i)(t/c)lV3 ♦ (toi.) 

To utilize the foregoing equations for actual computations, it is 
first necessary to evaluate Iw Since numerical values of ^ are 
available in the holograph only for -p = r| = 0, the evaluation must be 
carried out in two parts as follows: 


Iw = y n y ^ dTl' ( 4 l) 

-P « 00 

The first integral is evaluated from the results of the^numerical 

solution by mechanical integration of a curve of t) ''l'^('n,'^w) versus *t). 

6 

The values of the derivative used for this purpose are obtained from 
the equation 


e 



o 


( 42 ) 
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where the notation is as shown in the sketch. 

This equation is derived in the same way as equa- 
tion (3l)> except that is here taken equal 

to zero in accord with the boundary condition.^ 

It can be shown from Guderley’s singular solution 
for corner flow that for small negative values 

of T) the curve of rf must behave essen- 

Q 

tially as This result is useful in 

fairing the numerical results near t[ = 0. The 
first integral in equation (4l) contributes by 
far the majority (about 99 percent) of the total 
value of Iw • 

To evaluate the second integral in equation (4l) use in again made 
of the asymptotic solution (21). For this purpose, the constant B is 
determined such that the value of given by the asymptotic solution 

e 

matches the numerically determined value at the point ( -3>^w) • 
Substitution of equation (21) into the second integral of equation (4l) 
then gives 





exp 


■- -5r 

L 3 ^w J 


X 



exp 


jt 

3^w 


(-251)'^^=^ d?i 


(■‘S) 


where is 

0 

The integral on the 


determined from equation (42) applied at t] = -p. 
right is transformed through the substitution 


4- (-2?!)'’''= = “ 

3S« 


^he fact that the second derivative taken as zero in 

the present derivation follows from the boundary condition and the 
differential equation (l0)« 
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which gives 


r 
' 00 


exp 

jr 

3^ 

1 

/30W " 

v7/6 

00 


3X2®/^ V « / 


(-2ti) 3/21 d31 = 


u)^/6e-“du 


3©w 


The integral here is then found with sufficient accuracy by means of 
the following asymptotic formula, valid for large values of the lower 
limit (see pp. 95-96 of reference 3 I) : 



e”^ dll) = e”^ 


Equation (^3) thus reduces finally to 


With the value of Iv known, equations (37) and (4o) can be used 
to o'^ain the coordinates x/c ^and Y corresponding to any point in 
the t], 9 plane. The value of Y is obtained by direct substitution 
of the appropriate value of + into equation (40b). The value of x/c 
must be found by suitable integration of equation ( 37 ) • The location 
of the vertex A of the shock wave is found, for example, by integrat- 
ing equation (37) along the line OA in the holograph. (See sketch 
on p. 30.) If the leading edge 0 is taken as the origin in the 
physical plane, this gives 


M 00 

The integral here is evaluated in two parts following the procedure 
previously used in determining Iv 


For the abscissa of a point F on the shock wave, equation (37) 

gives 



(2^ d?i + de) 

e n 


where the integration is now taken from A to F along the shock 
polar. For purposes of numerical evaluation, the integrand here can be 
simplified by writing 


2 t ] dTi + \|f.^ d0 

e T) 


i|f 0 dff + i|r~ d© 
dTi + d9 

n 0 


2n+ (triA^)(d0/dr[) 
+ (d'^/d'ri) 

T) 0 


If and (d0/dT]’) are replaced by the appropriate functions of t( 

T) 0 

from equations (l^), there results finally 



(i^6) 


The inte gral in this equation is evaluated by plotting a curve of ^ 
versus yi+^ from the numerical results along the shock polar and 
carrying out the necessary integration by mechanical means. 


The abscissa of a point on the sonic line is found by integrating 
equation (37) along the 6 axis from B to C. Since point B is 
situated in the physical plane at x/c = l/2, this gives 



2 4lv T] 


(^7) 


The integral here^ is evaluated by mechanical integration of a curve of 
versus Q, where is found from equation (3l)* As can 

be seen from equation (28) and relation (29), tv.(O,0) in the vicinity 


of point B varies essentially as 
drawing the curve of ^^~(O,0) near 


. This fact is of use in 

/N/ ' 

0=0v It can further be seen with 


the aid of equation (40b) that near the shoulder of the wedge the trans- 
formed sonic line has the form 


\^Jc 2 

This relation is useful in establishing the detailed shape of the sonic 
line in the physical plane. It shows, in particular, that the sonic 
line will have a vertical tangent and an infinite curvature at the 
shoulder of the wedge. 
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Pressure Distribution and Drag 


To complete the analysis of the front wedge, it is left to deter- 
mine the pressure distribution and drag. Integration of eq^uation (37) 
gives for the chordwise location on the wedge of a given value of 


or 


5C 

c 


1 

2Iw 



0 





(48) 


The speed parameter | ^ (M^-l)/[(7+l)(t/c) which was used to 

present the results in reference 1, is related to rf by the following 
equation, derived from equations (20) and (21) of reference 1 and 
equations (l), (9)> (I 8 ) of the present paper: 


. 


2^-/3 
? 2/3 


(49) 


With these equations, the distribution of ^ as a function of x/c 
is readily determined. The integration of equation (48) is carried 
out by mechanical means using the same curve previously employed to 
determine Iw To fair the resulting 5 curve in the vicinity of the 
shoulder, use is again made of Guderley’ s analytical findings, which 
show that in this vicinity 



The chordwise distribution of the generalized pressure coefficient 







(r+i)^^°/(t/c)^^° 


Cp 


is related to the distribution of i by the following equation taken 
from pages l6 and 23 of reference 1; 


C p = - 2 (I-Iq) 


(50) 
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To find the portion of the generalized drag coefficient 
cd = [(7+l)/(t/c)^/®] cd 

contributed by the front half of the double-wedge profile^ equation (l3a) 
of reference 1 gives 


^df 




(51) 


To allow for the singularity in Cp at the leading edge, the integra- 
tion is best carried out directly in terms of the holograph variables. 
To this end, substitution from equations ( 31 ), (^9), and ( 50 ) gives 


~d. = 2 I, 


1 . L rf dq I 

iw j e J 


( 52 ) 


The integration here is performed in two parts following procedures 
parallel to those used in evaluating Iw 


CHAEACTERISTICS CONSTRUCTION OVER 
REAR OF AIRFOIL 

The characteristics in the plane (ti>0) are given by the 

fo].lowing relation obtained from equations ( 3 ) and ( 9 ): 

6 = const. ± rf 

3 


(53) 


The corresponding directions of the Mach lines in the generalized 
physical plane, as determined from this relation and the transformation 
equations (37) and (hOa),are 


dY 


d(x/c) 


= + 




(2n) 


1/2 


( 5 ^ 
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To the present order of approximation, therefore, the slope of the Mach 
lines is independent of the local inclination ?. This is consistent 
with the previous result that the streamlines must appear in the physi- 
cal plane as horizontal straight lines. As a consequence, the construc- 
tion of the Mach net over the rear of the airfoil is particularly 
simple in the small-disturbance theory. 

To in the construction, the equations for the characteristics 
in the t),0 plane can be conveniently written in the form 


9 = ^ 

3 




9 



rv» 'V* 

The symbols ©u and ©(^ denote, respectively, the ordinates at which 
the upgoing and downgoing characteristics through a given point (tj,0) 
intersect the vertical axis. Elimination of 6 between these equa- 
tions gives 




(55) 


which can be substituted into equation (5^) to obtain 


dY 


= 4 - 


/^ew/3\ V3 


± 


r\j 


■ 


d(x/c) \©(3-©uy 


(56) 


This is the basic relation for the characteristics construction in the 
physical plane. 

The construction of the Mach net itself follows a simple lattice- 
point procedure (cf. reference 32 )^ By identifying each Mach line with 
its appropriate value of ©u or the value of dY/d(x/c) at the 

intersection of any two Mach lines can easily be determined from 
equation (56) (or its graphical equivalent). The basic construction 
necessary to locate an unknown point c from the location of two known 
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points a and b is then as indicated in the sketch. The construction 
proceeds rapidly since, as pointed out, variations in the inclination of 



flow need not be considered in establishing the direction of the Mach 
lines. Where desired, the value of 0 can be found from the relation 



+ 0 


u 


2 


(57) 


The corresponding value of T] is given by equation (55).^^ 

Figure k shows a typical Mach net constructed by stepwise appli- 
£ation of the foregoing procedure. This net is for the case of 
^w (So = 0.921) and corresponds to the flow field shown for the 

front of the airfoil in figure 2 of reference 1. The construction is 


^*^In practice, the construction is actually carried out most easily in 
a plane of Y/(20^)^/^ versus x/c with the slope of the Mach 
lines given by 


d[~/(2a„)^/^] ^ ^ / 3/3 

d(x/e) 

This allows a single graph of slope versus to suffice for 

all values of ©w It also provides somewhat more convenient pro- 
portions for the construction of the Mach net. 
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begun at the shoulder of the airfoil (x/c = 1/2, Y = O) with the values 
of 0^ selected to provide approximately equal spacing between the Mach 
lines of the expansion fan. From the shoulder, the construction is 
carried outward to the sonic line and then inward to the rear surface 
of the airfoil. The drawing of the Mach line segments adjoining the 
sonic line might appear at first to offer some difficulty, since a 
linear average is obviously unreliable to determine a mean inclination 
in this vicinity. Actually, no trouble is encountered from this source, 
since the point at which each Mach line meets the sonic line is already 
known from the hodograph solution for the subsonic field. The construc- 
tion of the last segment approaching the sonic line thus reduces to a 
matter of simply connecting two known points. The slope of the first 
segment leaving the sonic line is found by either (a) multiplying the 
slope of the approaching segment by -1, or (b) determining a mean 
inclination based on the easily demonstrated fact that a Mach line in 
the vicinity of the sonic line behaves essentially as a semicubical 
parabola. It is immaterial to the final result which procedure is 
used. The identity of the Mach lines reflected from the rear surface 
of the airfoil is determined from equation (57) plus the boundary condi- 
tion that at this surface ^ As can be seen by comparing 

figure 4 of the present report with figure 2 of reference 1 , only a 
relatively small portion of the sonic line need be known to determine 
conditions on the rear of the airfoil. 


REMARKS ON ACCURACY OF SOLUTION 


Quantitative statements with regard to the accuracy of the present 
results are difficult to make. Fortunately, however, a check on the 
accuracy of the solution is available in the work itself. This check 
derives from the fact that, in the subsonic portion of the field, the 
calculated location of a given velocity in the physical plane should, 
theoretically, be independent of the path of integration which is 
followed in the hodograph. Thus, for example, the position of the 
velocity "t) = 0, = 1, which defines the point of intersection E of 

the shock wave and sonic line, should be the same irrespective of 
whether it is found from equations ( 45 ) and ( 46 ) 


(i) . j- 

\ C /e 2Iw 

or from equation (^7) 



e 



(58) 


11 


The latter possibility was pointed out to the authors by 
Gottfried Guderley. 
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t.(o,e) d0 

T1 


( 59 ) 


Actually, as was observed in connection with figure 1 of reference 1, 
the results of the two determinations show a small discrepancy. Such a 
result would be expected in any finite-difference solution. 


Discrepancies of the type noted can arise from two sources: (l) 

numerical inaccuracies in the relaxation solution of the finite- 
difference equations or in the transformation to the physical plane; (2) 
inaccuracies caused by the fact that the finite-difference equations 
themselves are not an exact representation of the boundary- value problem 
for the original partial differential equation. Experience with various 
refinements in the calculations indicates that the discrepancies here 
are primarily of the latter origin. Early computations with a coarse 
lattice and relatively crude finite-difference equations showed a con- 
siderable gap between the shock wave and the end of the sonic line. 
Increasing refinements in the grading of the lattice and in the deriva- 
tion of certain of the finite-difference equations gave progressive 
improvement in reducing this gap. This improvement came about primarily 
as a result of progressive reduction in the value of the integral 
the other integrals in equations ( 58 ) and (59) being relatively unaf- 
fected by the refinements in the calculations. Indications are that, 
in the results which were taken as final, the values of T>^w) 1 

G 

and hence of I^ are still somewhat too large. This means (see 
equation (40b)) that the ordinates of the shock wave and sonic line are 
probably somewhat smaller than they should be. The same is probably 
true, in general, of the corresponding values of |x/c|. Calculations 
of the chordwise distribution of ?) on the surface of the airfoil are, 
however, considerably more precise, since the errors in the two inte- 
grals in equation (48) tend to compensate. The refinements in the 
computations were, in fact, carried to the point where further better- 
ment caused only negligible change in the pressure distribution and 
over-all drag. Further evidence of the accuracy of the results in this 
regard is provided by the ease with which the computed values fair into 
the results of Guderley and Yoshihara at = 0 and into the analyti- 
cal curves which are available when the bow wave is attached and the 
flow is completely supersonic. (See figs. 4 and 6 of reference 1.) 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Oct. 8, 1951 
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TABLE I.- CALCULATED VALUES OF FOR = 1.6 Uq = 0.921) 
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Figure /. - Illustrative finite - difference lattice in the plane. 
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Figure 2. - Distribution of lattice spacing for ^ * L6. 
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